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ABSTRACT 


OIL RHEOLOGY ADJACENT TO THE SCRAPER RING 
OF A DIESEL ENGINE 


by 
DENNIS CHARLES LOGAN 


Submitted to the Department of Ocean Engineering in partial fulfillment of the 
requirements for the degrees of Master of Science in Naval Architecture/Marine 
Engineering and Master of Science in Mechanical Engineering. 


ABSTRACT 


Several recent experiments have been made to determine lubricant flow pattems 
in engine journal bearings. A laser fluorescence technique in use at the 
Massachusetts Institute of Technology allows accurate data collection of the oil film 
thickness on the ring pack of a production diesel engine. The data collected from the 
Kubota EA300N IDI engine consisted of five different types of lubricant--two single- 
grades, two multi-grades, and a synthetic multi-grade. 

The data was analyzed and it was found that while oil cross-flow 
circumferentially around the scraper ring is not present in fired cases, it is present in 
the motored cases. In addition, oil flow under the ring was evaluated with results 
consistent with previous observations and the flows observed suggest that the oil 
flow is between the ring and liner. Finally, a model for upstrokes that predicts the 
inlet wetting height for varying ring load was evaluated with positive results. 


Thesis Supervisor: Dr. David P. Hoult 
Title: Senior Research Associate 
Department of Mechanical Engineering 
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1.1 Background 

This research is the continued analysis of the lubricant film thickness in the 
piston ring region of a small production diesel engine. Olechowski [1] and Hoult, 
Wong, and Azzola [5] have recently investigated the various rheological effects of 
single and multi-grade lubricants in the top ring region. This project investigates the 
scraper ring (second ring) region using similar analysis techniques as used for the top 


ring analysis. 


One major goal of this research is to study fired data to expand the data base 
and to compare the results with the effects observed in the motored cases of the top 
ring. In addition, motored cases will be examined in the second ring analysis to note 
any anomalies between top ring and second ring observations. 

The data used for this analysis was collected using a laser fluorescence 
technique. McElwee [2] and Bliven [3] measured film thicknesses in the ring pack 
region using single-grade and multi-grade oils. McElwee collected data using five 
different types of oil -- two single-grades, two multi-grades and a synthetic multi- 
grade. Bliven's measurements consisted of one single-grade and one multi-grade of 
the same type that McElwee analyzed. In all cases, various operating conditions (ie. 
speed, load) were established to observe any variations in lubricant behavior as a 
function of operating conditions. 

The data was collected from a Kubota EA300N, a small, single cylinder 0.3 
liter indirect diesel engine. One major difference between McElwee's experiments 
and Bliven's experiments is the azimuthal location that oil film thickness was 


measured. McElwee measured oil thickness on the wrist pin axis while Bliven 





measured thickness on the antithrust side of the piston, approximately 65 degrees 
from McElwee's position. These two different measurements are important in the 
continued analysis of the oil flows in the piston ring region. The two different points 
of oil film measurement will play key roles in the determination of the boundary 
conditions upstream and downstream of each piston ring. 

The three phenomena investigated here all pertain to the scraper ring. They 


are cross-flow, flow under the ring and load versus inlet wetting height. 


1.2 Data Analysis Techniques 

The analysis of data begins with the ring fits on computer generated oil traces 
as shown in figure 1. The ring fitting technique is an art at best and is described by 
McElwee [2] and Olechowski [1]. After placement of the ring, measurements for 
numerical computations where taken from the fitted data. Figure 2 shows the various 
ring and film thickness parameters required for the calculations. 

One problem in analyzing the data was the fact that McElwee and Bliven used 
different film thickness calibration techniques [2,3]. This was accounted for by a 
temperature correction factor that was applied to McElwee's fired data. The 
correction factor was determined by a ratio of laser fluorescent efficiency at the 
various oil temperatures. The details of laser calibration are described by Hoult and 
Takiguchi [6]. McElwee's motored data and all of Bliven's data did not require a 
correction due to the dynamic (self-calibrating) calibration used in these cases. 

The various constants and oil parameters used were taken from previous 
work and actual measurement. The scraper ring geometry was established using 
talysurf measurements of the actual ring used in both experiments. Appendix A 
describes the talysurf procedure used in determining the scraper ring parameters. 

The oil temperatures used in the fired analysis were obtained from Hartman [4] who 


determined cylinder liner temperatures under loaded conditions. With an oil 





thickness of only a few microns, it is assumed that cylinder liner temperature is the 
most accurate measurement of oil temperature in the ring pack region. Motored oil 
temperatures were the same used by Olechowski [1]. 


The data base analyzed is shown in appendix B. 
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CHAPTER 2 - THEORY 


2.1 Cross-flow 


Cross-flow is the most straight forward of the three phenomena to analyze. 


A comparison of the minimum film thickness ho between the two azimuthal positions 


was used to determine the presence of cross-flow. Since approximate calculations 


(appendix C) showed that ho/h.. ~ 2, the comparison of hy between positions is a 
proportional comparison of hs at each location and thus gives an indication of flow 


from one position to the other. 


2.2 Flow Under the Ring 
Flow under the ring begins with the Reynolds flow equation as defined by 
Coyne and Elrod [7]: 


The boundary conditions here are: 


1. ring load (W) 
2. film thickness "far" away from the ring (hoo) 
3. ho or separation point at the exit (X2) 


4. known ring profile. 


The nondimensional flow equation is obtained by using the standard 


nondimensional terms h, x and P: 





h= 





P= Ph? 
Uub 


på 
li 


ehe 
ho 


০1৯4 


By substitution, the flow in nondimensional terms becomes: 





From conservation of mass (no sources), dQ/dx = O or Q = a constant. Also, 
because of the relatively high velocity of the ring relative to cross-flow velocities, it 
seems a very good approximation to set Q = 1 under the ring. Therefore, the 


relationship between the entrance point (x =1 ) and exit point (x = 0) is 


h° 


I 
A 
= 


(1 P(1) = h(0) এ 


টু. 12 NOU P(0) 


h(0) _ 6- R0) PO) 
hl) 6-10) P(0) 


and the relationship between the entrance film thickness and hos is 


KD h (1) dP(1) 
12 dx. 


Flows under the ring were analyzed using the same computer programs as 


Bliven [3] and Hartman [4] used in their analysis. 


m N 





2.3 Load versus Inlet Wetting 

Similar to the top ring analysis [1,5], the variation in load versus inlet wetting 
height was investigated. The nondimensional load G is the same used by HWA [5] 
but normalized for the scraper ring (ie APB for the scraper ring). G, the bearing 


number, is defined as follows: 


2 
G = Ub 
h2APB 





₹ 0071, 12, P2, P3) 


where P2, P3 are the nondimensional pressures of the second and third lands 
respectively. To eliminate some of the variables it was determined that the pressure 
differences across the ring were small enough to disregard P2 and P3 for the 
compression and exhaust strokes. In addition, the assumption was made that I? 
always occurred, on upstrokes, at the sharp undercut tip of the scraper ring [9] with 
ho and I determined based on inputs b, I2. The algorithm used to determine G for 


upstrokes is described in appendix C. 
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3.1 Cross-flow 

Figures 3 and 4 show the plots of the wrist pin position ho (McElwee) versus 
the skirt position ho (Bliven) for a typical motored and fired case respectively. 

Figure 3 shows that the average value of ho (with standard deviation) is mostly above 
the 1:1 line which indicates that cross-flows are probable in the motored cases. 
Figure 4 shows a similar plot for a fired case with the values much closer to the 1:1 
line. This is an indication that cross-flow is unlikely in the fired cases. 

These results help define the boundary conditions in that differences in 
observed and calculated h.. could be due to cross-flows around the piston. 
Preliminary observations, in fact, show a difference of a factor of approximately two 
between h.o observed and h.o calculated. The cross-flow results provide more clues 
to the correct boundary conditions required to solve the flow problem. Additional 


plots are shown in appendix D. 


3.2 Flow Under the Ring 

Figure 5 shows the difference in oil flow rates between power 
(compression/expansion) and gas exchange (intake/exhaust) strokes for four fired 
and four motored cases along with the previous observations [3,8]. The results are 
consistent with previous observations which conclude that distinct oil recirculation 
loops appear to exist in the piston ring region and between the piston and the sump. 
The flow of oil under the second ring is part of the smaller ring region loop (the 
major loop is through the oil relief holes under the oil control ring that direct oil back 
to the sump). In addition, the oil transport rates are on the order of 10 to 100 times 


larger than oil consumption rates [4] which is consistent with previous observations. 
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There are three paths by which oil can reach the second land: flow through 
the ring gaps, flow around the backside of the ring and flow between the ring and the 
liner. Only flow between the ring and the liner is observed in this work. The fact 
that the flow between the ring and the liner, presented here, agrees with the flows to 
(from) the second land [8] strongly suggests that for the second land, the lubricant 
flow path is between the 2nd ring and the liner. 

Additional condensed plots of flow segregating power strokes, gas exchange 


strokes and overall flow for the second and third lands are shown in appendix D. 


3.3 Load versus Inlet Wetting Height 


Figure 6 shows the correlation of the theoretical upstroke model with 
experimental data for fired cases. The model exhibits two branches which 
approximates the observed data fairly well. Figure 7 shows the plot of I'; versus I? 
as a correlation with the upstroke model. The values of I > greater than one agree 
with the model to a high degree while the values close to one do not agree. For the 
values close to one, the outlet separation point is moving away from the undercut tip 
toward xo (h(xo) = ho) and is behaving similar to the top ring's outlet separation 
point. This is because above the undercut, the scraper ring has a circular profile 
similar to the top ring and thus top ring theory will apply when I> is close to one. 


Similar plots for the motored case are shown in appendix D. 
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The results of chapter 3 lead to the following conclusions concerning oil flow 


adjacent to the scraper ring: 


1. Little or no cross-flow present in the fired case and probable cross-flows present 
in the motored case. These results provide a better understanding of the observed 
film heights observed away from the ring which allow more appropriate boundary 


conditions to be defined. 


2. Oil flows under the ring are part of the ring region oil transport loop and are much 
larger than oil consumption rates which is consistent with previous observations. 
Also, the present results suggest that the flow to (from) the second land arises from 


flow between the ring and liner. 


3. The relationship between load and inlet wetting height, as determined by the 


proposed upstroke model, is in good agreement with lubrication theory [9]. 
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Figure 1 - General scraper ring fit on a computer generated oil trace 








Lubricant hee ~ ho/2 
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E a | a b 2 | Liner 


Figure 2 - Definition of terms used in the solution of Reynolds flow equation under 


the scraper ring. 
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Figure 3 - Typical plot of minimum film thickness (ho) at each azimuthal location for 
a motored case (pen multi 3000 rpm). The average values indicate probable cross- 
flow from the skirt position to the wrist pin position in this motored case. 
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Figure 4 - Typical Plot of minimum film thickness (ho) at each azimuthal position for 
a fired case (pen single 1500 rpm). The average values indicate unlikely cross-flow 
in the fired case. 
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Figure 5 - Typical plots of oil flow rate comparisons of flow under the second ring 
for power (PWR) and gas exchange (GAS) strokes in fired and motored cases. 
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Figure 6 - Plot of nondimensional load G versus I; for fired upstrokes. The solid 


line represents the theoretical upstroke model which assumes the exit separation point 
(X2) is located at the undercut tip of the ring. 


21 








FIRED UPSTROKES 


hav single 
pen single 
pen multi 
cum multi 
mob multi 
— — Theory 





Figure 7 - Inlet height versus outlet height with the line representing the theoretical 
upstroke model. 
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Appendix A - Scraper Ring Talysurf Profile/Ring Tip Radius 
Calculation 
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RING TIP RADIUS CALCULATION 


The scraper ring tip, above the undercut, can be fit with a circular profile using a 


Taylor Series expansion of the equation of a circle around ho which gives: 


_ (x-xo)? 
h(x) = ho m (1) 


where a is the radius of the circle. Since the arc of the circle used to describe the ring tip 
is small, it can also be described by the tip of a parabola. Therefore, expanding equation 
(1) and setting it equal to the oil film height equation as a function of the the ring tip 
(from the Talysurf plot), a can be found by matching any same term of each side of the 


equation. An example using the talysurf on the previous page is shown below. 









h=ho+(s-y) 
s=maxy 
= 22.8 um 


y = 22.6 + 2.4x - 7.7x42 


xo = .153 mm X 
y —CI + ০2% + C3X2 


(x2 + 2xox + ৯0) 


h = họ + (22.8 - cı - ০2% - c3x?) =ho + 25 


Matching the squared terms, we get 


-C3 — —- 


or, using the parabolic fit terms, a — .065 m. This is the ring tip radius of a circular tip 
profile of the scraper ring. 
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Appendix B - Data Base Analyzed 
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The data base used in this research was generated from experimental data taken 
by McElwee [2] and Bliven [3]. The data set numbers designate, among other things, 
which person took the data. Data sets numbered less than 500 were taken by McElwee 
and sets numbered 500 or greater were taken by Bliven. The cross-flow analysis data 
base used two oils at two different speeds each since Bliven only worked with two oils 
which were a subset of the five oils analyzed. The load versus inlet wetting height data 
base was generated using a cross-section of all five oils at five speeds. The cross-flow 
data base is shown in table 1. The load versus inlet wetting height data base is shown in 
table 2. 

The actual data base is shown in the rest of appendix B. Each data set is stored in 
an Excel spreadsheet and saved on floppy disks. The name used to save each data set is 
coded as follows: 


Character # Translation 
] ring (S=second) 
2 data (B=Bliven, M=McElwee) 
3 oil (P=Pennzoil, H=Havoline, C=Cummins, M=Mobill) 
4 type (S=single-grade, M=multi-grade) 
5,6 speed (add two zeros to this number, ie 25=2500) 
7 load (H=half load, F=full load). 


Motored data sets are also distinguishable by a ".X" at the end of the name. 

Two different types of measurements were used when collecting the data. 
Initially, a digitizer was used to measure each point on a particular oil trace. This is 
shown by columns in a particular spreadsheet labeled x1, x2, x3, yl, y2, y3 (1=inlet, 
2=minimum oil height points, 3=exit). The x values (distance along piston) are mm and 
the y values (oil height) are um. These values were taken straight off of the oil traces and 
entered as dimensional values. Since this procedure required extensive amounts of time, 
a second method was devised using an engineer's ruler on the 30 divisions/inch scale 
measuring the number of tick marks (divisions) for a particular point and then scaling the 
length with the x or y axis as given on each oil trace. The spreadsheet columns using 
this procedure are designated b ticks (only one measurement required to obtain b), hl, h2 
or h3 ticks (the oil height measured in ruler divisions where 1,2,3 have the same 
meaning as above). These values were then scaled according to the x or y axis scale. 

The various fired and motored data constants used in the calculations are stored in 
separate spreadsheets. These are shown at the end of appendix B. 
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Fired (all full load) 


TEGE BE TT ETE 


Motored 


Sa S USUS NS 


Data Set # 


359 
368 
409 
418 
521 
522 
526 
527 


392 
355 
402 
405 
519 
520 
524 
525 


Oil 
Pennzoil 15W40 
Pennzoil 15W40 
Pennzoil SAE30 
Pennzoil SAE30 
Pennzoil 15W40 
Pennzoil 15W40 
Pennzoil SAE30 
Pennzoil SAE30 


Pennzoil 15W40 
Pennzoil 15W40 
Pennzoil SAE30 
Pennzoil SAE30 
Pennzoil 15W40 
Pennzoil 15W40 
Pennzoil SAE30 
Pennzoil SAE30 


Table 1 - Cross-flow Data Base 
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1500 
3000 
1500 
3000 
1500 
3000 
1500 
3000 


1500 
3000 
1500 
3000 
1500 
3000 
1500 
3000 





Fired 


ST]! A A ই 


Motored 


SET s 


Dat 
207 
258 
263 
268 
306 
pal 
522 
526 
527 


203 
Zl 
254 
303 
519 
520 
524 
525 


Table 2 - Load versus Inlet Wetting Height Data Base 


Oil 


Cummins 15W40 


Havoline SAE30 
Havoline SAE30 
Havoline SAE30 
Mobill 15W50 

Pennzoil 15W40 
Pennzoil 15W40 
Pennzoil SAE30 
Pennzoil SAE30 


Cummins 15W40 


Havoline SAE30 
Havoline SAE30 
Mobill 15W50 

Pennzoil 15W40 
Pennzoil 15W40 
Pennzoil SAE30 
Pennzoil SAE30 
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Speed (rpm) 
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1000 
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1500 
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1500 
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2000 
1500 
3000 
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Appendix C - Upstroke Model Development 
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Definition of terms used for profile of second ring as defined by Coyne and Elrod [7]: 


[1 - nondimensional inlet wetting height 
I5 - nondimensional outlet wetting height 
ho - minimum wetting height 

X1 - local position of I y relative to 0 

O - local position of ho 

X2 - local position of I^ relative to 0 

a - radius of parabolic fit of ring tip 

b - wetted length 

05 - I? - ho 


The formula describing the oil height under the ring is 


h(x) = no + ZZO 
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where xo = O (local position of ho). Based on inputs b, I? with a, 62 and X2 known 
from the scraper ring talysurf calculations, ho and I | can be determined using the above 
equation. 

These parameters are used in an Excel macro called "Macro 9". Macro 9 is an 
iterative scheme to compute f, [ 1, b and h.. from a, X2, I 2 and W (load). X2 is the 
separation location measured from the point on the ring where h = ho. There is a main 
program and five subroutines called K, Ksubl, hinfh0, loc1 and h0. The subroutine 
hinfh0 calculates hinf/ho. K and Ksubl calculate the values of K in the formulae of 
Coyne and Elrod while locl calculates the value of XI (the distance under the ring from 
the inlet to the point where h = ho, nondimensionalized). The subroutine h0 calculates 
ho. Main calls hO which calculates ho directly from I», X2 and a. Then the 
nondimensional variables for X2 and W are calculated. The accuracy constant, called 
eps for epsilon in the code, is assumed at 0.0001. The iteration variable is X1 which is 
looped in main through subroutine locl. Locl calls three subroutines, hinfh0, K and 
Ksubl. Each uses the formulae given in Coyne and Elrod [7]. Locl is a root solver 
using a bisection technique to guess the most accurate X1, the root of the equation for 
Wcalc - W, where Wcalc is the load calculated from all other variables. The bisection 
technique progressively narrows the location of the root down until the accuracy criterion 
is met. The last guess of X1 is then sent back to main where the entire loop is done 
again until X1 converges to the accuracy limit. 

Macro 9 was written by Jim H. Azzola, PhD candidate, MIT 1991, Sloan Auto 
Lab. 
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MPs hos — ee | 
MEROS SS A es 
Eee :5- রা HCH I 
| 3 |mein 

| 4 |-RESULT(64) 

| 5 |[=ARGUMENT("mu", 1 gam2' 

6 |=ARGUMENT('W) 000 00 0 0 


| 8 |-ARGUMENTCs,D 0. | বৰ =xr"2/(2*a*(gam2-1)) 
sARGUMENT("xr",1) 


|11 JeSET.VALUE(A12,MACRO9. XLM !hO(gam2,,x1)) ST PT 
| 12 |0.000017777777777778 O a) 
QRT(2*a*A12))) (set non-dimensional x2 | sid 


set non-dimensional x2 












main program indirect calc | 


























| 15 |-SET. VALUE(A16,"w*A12/(12*mu*u*a)) 


| 16|0.081705948372615 ^ [w*snondimensionalloed |, || 
PRB | rt [| |) | 
| 18 |- SET. VALUE(A20,-0.7) setinitialxlguess | | 
| 19 | ET XEM 
|20 |-0.8999267578125  frew | | | | O 
de —_ — — — XM ছি! 
TAI A o O ee a a 
| 24 | FSELVALUBAMOMACRONXLMISENAISAISATD) pet new valde U — I 
EXER I] ma * ।৷। 

28 1-0.8999267578125 —  — fold 1 | 
LE ASSES big হী 
=IF(ABS(A20-A26)<=A 17,GOTO(A29),GOTO(A22)) ES pc 0 TE 
| 30 |=A29*TAN(A20) Dimensional left value x1 







OTER ESSEN a erg 
FS2[-SETVALUEQGRAI) (mw 1 
[33 |-SET. VALUE(B38,A30 সি E 
734 |-SET.VALUE(C38 MACROS XLMIBRBOAZO AT AINE  |int 1 TT 
8 j-RETURN(A3S:C4) rearme vas T_T] 
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[4 |-RESULT() 
|. 5 |ZARGUMENT("exr^,1) 


; 


initial right value of interval 
=SET.VALUE(D22,D10) ভি X 
=SET.VALUE(D23,D11 Ber 


| 14 I=SET.VALUE(D24,0) 
| 15 |=(D22+D23)/2 | 
| 16 |=IF(D23-D 15<eps,GOTO(D27),GOTO(D17)) 


zero counter 
ND SE 
17 |=-SET.VALUE(D18.MACRO9.XLM !hinfhO(D15,exr ভেড়া 


1.1662557645283 *hinf/h0 = T from C&E 
| 19 |[=MACRO9.XLMIK(D15,D18)-MACRO9.XLMIK (exr,D18+MACRO9.XLM!Ksub1(D15,D18)*(TAN(exr)-TAN(D15))-wstar|f(x left int val) 


MACRO9.XLMIK (D23,D18)-MACRO9.XLMIK (exr,D 18 +MACRO9.XLMIK sub 1(D23,D18)*(TAN(exr)-TAN(D23))-wstar f(x right int val) i 
D20)<=0, SET. VALUE(D22,D 15), SET. VALUE(D23,D 15) 


left interval val 


right int value 
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| 


_____|Ksubl 

=RESULT(1) =RESULT(1) 

= mARGUMENT("ex",1) 
=ARGUMENT("ef",1 =ARGUMENT("ef",1) 
=-ex*TAN(ex)/2+ef*(3*ex* TAN(ex)/8-COS(2*ex)/16) |=-ex/2-SIN(2*ex)/4+ef*(3*ex/8+SIN(2*ex)/4+SIN(4*ex)/32 
=RETURN(F7 zRETURN(G? 
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[4 |=RESULT(1) mu 

| 5 |=ARGUMENT( exl",1) = 

=ARGUMENT("exr" 1) 
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=ARGUMENT("b", I 
=ARGUMENT("hsubO", 1 





| 10 |ZARGUMENTN(" hint", m 
(11 JSARGUMENT(”sv”,D — 

(12 JSARGUMENT('xl”,) [dim zi 
| 13 |mARGUMENT("xr",1) dim x2 
(14 |-SQRT(2*a*hsubo) 

| 15 |=ATAN(x1/114) 

17 |=I15-(0.5*TAN(15){1+TAN(15))}X2}+I15)*0.5*hinf/hsub0 
=116-{0.5*TAN(116)/(1+(TAN(116))A2}I16)*0.5*hinf/hsub0 

| 19 [=mu*u*4*a/114*(118-117) 

| 20 {=I19/(u*mu*b/hsub0 
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Appendix D - Additional Figures 
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LIST OF FIGURES 


Figures 3A, 3B, 3C: Cross-flow plots - motored 
Figures 4A, 4B, 4€: Cross-flow plots - fired 
Figures SA-G: Flow under the ring plots 

Figure 6A: Upstroke model G versus I] - motored 


Figure 7A: Upstroke model I ] versus [2 - motored 


KEY 
pm - pen multi 
ps - pen single 
15 - 1500rpm 
30 - 3000rpm 
M - Wrist pin (McElwee) 
B - Skirt (Bliven) 
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